Quantum computation with doped silicon cavities 
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We propose a quantum computer architecture involving substitutional donors in photonic-crystal silicon cav- 
ities and the optical initialization, manipulation, and detection processes already demonstrated in ion traps and 
other atomic systems. Our scheme considerably simplifies the implementation of the building blocks for the 
successful operation of silicon-based solid-state quantum computers, including positioning of the donors, real- 
ization of one- and two-qubit gates, initialization and readout of the qubits. Detailed consideration of the pro- 
cesses involved, using state-of-the-art values for the relevant parameters, indicates that this architecture might 
lead to errors per gate compatible with scalable quantum computation. 

PACS numbers: 03.67.Lx, 71.55.Cn, 42.50.Pq 



The search for a working quantum computer has comprised 
areas ranging from optics to atomic and condensed-matter 
physics [I J. Finding physical systems that allow for accurate 
operations has been a formidable challenge, yet to be met. 
Indeed, the viability of quantum computers depends on find- 
ing physical systems that allow scalable fault-tolerant com- 
putation, which means that the errors remain bounded when 
the number of qubits increases. In order to have scalable 
quantum computation, the error per gate (EPG) should be 
smaller than a certain threshold, which depends on the spe- 
cific error-correction scheme. For independent and identically 
distributed errors, the best lower bound so far is 1.9 x 10""* 
(12 • For other architectures, which require however a large 
resource overhead lH], this threshold is bounded below by 
1.04 X 10"^ iH. 

Semiconductor devices ||5]|6l, and most particularly those 
based on silicon [7|, have attracted considerable attention, but 
actual realization is hindered by difficulties concerning scal- 
ability, detection and fabrication fS). Most candidates for 
a semiconductor-based quantum computer rely on spin- 1/2 
fermion qubits |9|, which for Si may be associated to the 
long-lived electron and nuclear spins of shallow substitutional 
donors |7|. In fact, electronic spin decohere nee times have 
been shown to be larger than 60 ms in ^^Si (isotopically pu- 
rified) at temperatures of 7 K |10|. Implementation of quan- 
tum computation with these systems is hindered by several 
problems. The most obvious ID is the difficulty in the ma- 
nipulation and measurement of single-spin states. Two-qubit 
operations relying on exchange gates |5|, restricted to nearest- 
neighbor interactions, have limited scalability potential ifTTl . 
The particular electronic band structure of bulk Si leads to fast 
oscillations in the electronic exchange coupling when the in- 
teracting donor-pair relative position is changed on a lattice- 
parameter scale |12|. Thus, proposals based on this mech- 
anism require nanofabrication techniques far beyond current 
capabilities. 

Here we propose placing the donors in a single-mode 
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FIG. 1 : Proposed quantum-computer architecture. Donor impuri- 
ties are placed in the neighborhood of intensity maxima of a photonic 
crystal cavity mode, not necessarily every maximum. The donors are 
under the action of a uniform magnetic field B and electric fields E, 
produced by the electrodes E. The magnetic field is strong enough to 
decouple nuclear and electronic spins in the donor ground state, and 
Zeeman-splits two electronic spin states, j and | [), which con- 
stitute the qubit. Turning on the electric fields allows to switch off 
individual qubit Raman transitions induced by the two laser beams, 
spread out over the ensemble of qubits, and also the coupling among 
different qubits through the vacuum cavity field. The inset displays 
the misplacement 5r of an impurity from a maximum of the cavity 
mode. 



photonic-crystal Si cavity flSl, and optically addressing them 
through Raman transitions induced by the cavity mode and 
only three applied laser beams, spread out over the whole en- 
semble. A fourth laser beam is used for the readout. The 
qubits result from the interaction of the electron spin with a 
uniform magnetic field {E). The system is kept at a tem- 
perature around 7 K. Essential elements of the architecture 
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FIG. 2: Orbital energy levels of As in Si. Only the relevant levels 
are shown. Here CB stands for conduction band. 



considered here are schematically illustrated in Fig. [T| where 
we represent an array of donors positioned at the maxima of 
the cavity mode. One- and two-qubit logical gates, as well 
as system initialization and readout, are implemented through 
the external laser beams, which address all qubits simultane- 
ously. The coupling between a donor electron and the light 
fields may be interrupted by an external electric field, due to 
the Stark shift of the donor levels. This effect is explored for 
selecting the target qubits for one- or two-qubit operations. 
Two-qubit operations are mediated by the vacuum field of the 
photonic-crystal cavity |6|. 

The study of group-V donor impurities in silicon is a quite 
mature field llT4l [TSll . When a group-V element, such as P, 
As or Sb, substitutes the group-IV Si atom in bulk Si, an ad- 
ditional electron is incorporated in the system. This electron 
remains bound to the core potential via a screened Coulomb 
interaction, constituting a solid-state analogue of the hydrogen 
atom. Electronic-structure peculiarities of bulk silicon lead to 
a modified hydogenic spectrum for the donor. Experimental 
values for the lowest energy levels and relative energies for As 
donors in silicon [15] are given in Fig. [2] 

Qubit states are defined within the 15(^41) (see Fig. |2]i 
ground state under a magnetic field B satisfying gefisB ^ 
A, where A is the hyperfine coupling constant, fiB is the Bohr 
magneton and .ge(~ 2) is the electron Lande factor in Si. For 
As in Si, A = 400 MHz, so a magnetic field B = 0.3 T will 
decouple nuclear and electron spin states, generating well- 
defined electron spin states | J,) and | f) as qubits, with a 
10 GHz splitting. 

One-qubit operations are implemented through Raman cou- 
pling of the states | I) and | f) of donors previously se- 
lected by switching off the Stark-shift electric fields acting 
on them. In this excitation scheme, the donors interact with 
two laser fields of Rabi frequencies J^LiI and ^Il2^ and fre- 
quencies and ojl-,, respectively, detuned from the tran- 
sitions between the qubit states | j) and | f) and the states 
|2P+__) by A, as shown in Fig.[3|a) (the fine structure split- 
ting of the levels 2P+ _ is not shown). If A ^ fi^j ,^L2i Tp, 
with Tp being the decay rate of the levels 2P_|_ the lev- 
els 2P+ _ are only virtually populated, giving rise to an ef- 



fective coupling between levels | j.) and | "f), described, in 
the interaction picture, if = i^i2' by the Hamiltonian 
i/eff = ^^ff U)(T I + H.C., where = jf^ij' /A. 

During the Raman coupling of the qubit states there is a 
small probability, of the order f^off / A, of populating the inter- 
mediate level 2P+ _ . This will lead to decoherence of the one- 
qubit operations with the rate fioff (Fp/A) due to the sponta- 
neous decay of the level 2P_|_ Since the time required for 
one-qubit operations is of the order of l/ficff, the error prob- 
ability per gate will be ei w Tp/A. 

We propose a spin-orbit (SO) mediated coupling for the 
opposite-spin qubit states. In order to avoid destructive in- 
terference effects that would make ficfi vanishingly small, the 
intensity C, of the SO coupling among the states within the 
2P_|_ _ manifold must be comparable or larger than the de- 
tuning A of the laser fields. Strong SO splittings have been 
measured for states of the fundamental manifold in the group 
VI donors Se:Si and Te:Si, and of their corresponding ionized 
states ifT^ [TTl [181 . For Se:Si, the measured SO coupling is 
3.2 cm~^ (96 GHz); so we assume in our calculations for As, 
the element corresponding to the Se row in the periodic table, 
a coupling of the same order of magnitude, ^ 100 GHz. 

The SO coupUng also mediates an efficient mechanism 
for selectively turning on and off the interaction between the 
qubits and the light fields through an applied electric field pro- 
duced by the electrode above each donor. The electric field 
has two effects: It increases the detuning between the Raman 
laser fields and the atomic transition, so that it becomes much 
higher than the SO splitting, and it mixes 2P and 2S states. 
The increase of the detuning causes a destructive interference 
between the 2P states, which leads to the vanishing of the Ra- 
man transition. The mixing of 2P and 2S states also reduces 
the Raman coupling, since the 2S state does not couple with 
the ground state. Perturbation theory indicates that the com- 
bined effect reduces the transition probability by a factor equal 
to the third power of the ratio between the SO coupling and 
the electric dipole energy, which leads to an error of the or- 
der of 10~^ for an applied field equal to 20 kV/cm. This field 
is below the ionization threshold for P in Si 1 19 1, and for As 
the threshold shoud be even higher since the binding ener- 
gies are larger Only the donors that are subjected to smaller 
electric fields will be affected by the Raman coupling. This 
scheme has the advantage that only quiescent atoms are sub- 
ject to large electric fields, so that the essential properties of 
the active atoms remain unchanged. 

One-qubit operations are implemented with a linearly- 
polarized beam along the direction of the magnetic field (same 
polarization as for the cavity mode) and a circularly right- 
polarized beam, as shown in Fig.[3]^a). Due to selection rules, 
the cavity mode does not affect one-qubit operations. 

Two-qubit operations involve the interaction of a previously 
Stark-shift-selected pair of qubits with an additional laser 
beam and with the cavity mode. The frequencies and polar- 
izations of the laser and cavity fields are chosen in such a way 
that they nearly satisfy the conditions for Raman coupling of 
the qubit states | I) and | f) for each donor |6| (see Fig.jsjb)). 
The wavelength for the transition lS{Ai) 2P_|_ _ is 26.39 
fim in vacuum; in Si (dielectric constant e = 11.4), the cor- 
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responding value is A = 7.82 /im. The wavelengths of the 
cavity mode and the laser beams should be around this value. 

Under the conditions established for a Raman transition and 
in the dispersive regime, i.e., ^ ilL,Tp and ^ 
r^c, Tp, Fc, where ilc quantifies the coupling with the cavity 
mode, with width Tc, the states 2P+__ are only virtually oc- 
cupied, giving rise to an effective coupling between the cavity 
mode and the qubit states which, in an adequate interaction 
picture, is described by the Hamiltonian: 
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where a is the annihilation operator for the cavity field, — 
I i)ii(t I is the spin flip operator for donor i, 6^ — A\ — A^, 
and the sum extends over all the donors selected by the electric 
static fields. The couplings fi^g are defined as: 
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If 5i ^ ri*ff , Fc, the cavity field will be only virtually ex- 
cited and can be eliminated from the dynamics, leading to an 
effective two-qubit interaction mediated by the vacuum of the 
cavity mode: 
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H.c. 



(3) 



where J'j = (5* — . From Eq. one can see that each 
pair of qubits i, j that satisfies 5* = 5^ will resonantly interact 
through the Hamiltonian 



with 



the effective 



H.c. 



constant 



(4) 



will 



coupling ai,ij 

/(5*. The qubit pairs for which S^' ^ fiij 

interact off-resonantly and will not couple to each other. The 
error probability per gate for two-qubit operations can be 
found in a similar way as for one-qubit operations and will be 
given by £2 ~ T^c/S'. 
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FIG. 3: Logical gates, (a) One-qubit operations are implemented 
with a linearly polarized beam along the direction of the magnetic 
field and a circularly right-polarized beam; (b) Two-qubit operations 
are implemented with a linearly polarized cavity field along the direc- 
tion of the magnetic field and a circularly left-polarized laser beam. 




FIG. 4: Initialization and readout, (a) Initialization: population of 
qubit state | t) is transferred to the qubit state j I) via optical pump- 
ing by resonantly exciting transition | |) <-> 2P+,_ with laser light; 
(b)Readout of the qubit-state is made by monitoring the fluorescence 
light of the cyclic transition | T) ^ 2Po. 



The Hamiltonian Q implements the VSWAP operation 
I Ti) (I U) + I iT))/%/2, which, combined with single- 
qubit rotations, can be used to implement a CNOT gate |5|. 
VSWAP operations can be implemented in parallel, by hav- 



ing different pairs. 

I a,- 1. 



with ^ S^, = S\ but \S' 



> 



The qubits are initialized by driving resonantly the transi- 
tion I t) ^ 2P+ _ with a laser beam. Under the action of the 
magnetic field B and at temperatures of the order of 7 K only 
the state 1S{A1) will be populated. Since the level 2P_|. „ is 
unstable, it will decay to one of the qubit levels, leading to 
optical pumping of the level | |) (see Fig.|4];a)). 

Qubit readout takes advantage of the fact that the states of 
the 2Po manifold do not show SO coupling. If laser light ex- 
cites resonantly the transition | f) ^ 2Po, only the states of 
that manifold with the same electronic spin as the state | t) 
are excited. Due to selection rules, the radiative- or phonon- 
assisted decay of these states to states with different electronic 
spin is forbidden. The decay out of the 2Po level is radiative, 
whereas the decay out of level IS{E) is phonon assisted. For 
this reason, all the excitation will decay back to the state | f). 
Therefore the transition | f) ^ 2Po is cyclic and the elec- 
tron shelving technique can be used to measure the occupa- 
tion of the qubit states |20|: If, during the laser excitation of 
the transition | f) ^ 2Po fluorescence light is observed, the 
qubit was in state | f), otherwise the state | |) was occupied 
(see Fig. |4|b)). Since the decay 1S{E) ^ | t) is assisted 
by acoustic phonons, part of the fluorescence light differs in 
frequency from the laser exciting the transition | "f) ^ 2Po, 
which implies that it is possible to distinguish the fluorescence 
signal from scattered laser radiation. 

The feasibility of the proposed scheme is based on the fol- 
lowing estimates for the frequencies, couplings, and times in- 
volved in the one-qubit, two-qubit and readout operations. 

The measured absorption linewidth of the IS(Al) — > 
2P^ _ transition for Si:P is approximately IGHz |21 1, giv- 
ing an upper bound for the decay rate Fp. This rate could be 
strongly decreased (more than one order of magnitude), since 
phonon-mediated decay can be suppressed by applying stress, 
as demonstrated in Ref. |22|, and spontaneous radiative tran- 
sitions from these levels are also strongly supressed due to the 
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photonic band gap, since they are far detuned from the cavity 
mode of the photonic crystal [23 1. The Raman-coupling con- 
ditions are satisfied, for example, for riij/27r = ^1^^/211 = 
2GHz and A = 200GHz. This would lead to an effective 
Rabi frequency ricff/27r — 20MHz and an error probability 
per gate ei w 10^**. For these parameters, which correspond 
to laser powers of the order of 10 mW, the typical time for 
a one-qubit operation would be of the order of 50 ns, much 
shorter than a spin decoherence time of 60 ms. Under the 
same conditions, we have calculated that the error probabil- 
ity per gate induced by eventual impurity ionization, due to 
two-photon absorption, is neglig ibly small (e w 6 x 10""^). 

For the two-qubit operations, one could choose, for exam- 
ple, Q.\j/2ti ^ 30MHz (this value corresponds to a modal 
volume lOOA^), n'^/2TT - 5GHz, A*^ = lOOGHz and 
A^ = 99GHZ. This yields an effective two-qubit coupling 
f2ij/27r ~ 2.25KHz, which allows one to perform more than 
10^ VSWAP operations within a qubit decoherence time of 
60 ms. Since for two-qubit operations the error per gate is 
£2 ~ Tc /Si, an error per gate of the order of 1 x lO^'^ would 
imply a cavity decay rate Vc ~ IMHz. This requires a cavity 
quality factor Q k, 10^. Quality factors of 10® have already 
been reported for silicon-based photonic-crystal nanocavities 
|U3l; Q's as high as 2 x 10*^ seem to be within reach |13|. 
Larger wavelengths in the fim region, as used in our proposal, 
should lead to yet higher values of Q. Combined with larger 
values of the spin decoherence time, consistent with the ex- 
perimental results ifTOl l24l . this would allow one to increase 
5, further reducing the error per gate. 

Readout is very fast, since the decay rate out of level Pq is 
of the order of 1 GHz whereas the phonon-assisted decay of 
level IS{E) is of the order of 10^^° s ||25l. We estimate that 
about 10000 cycles, corresponding to a detection time around 
1 /xs, should yield a reading efficiency close to 100%. Parallel 
readout can be implemented for donors separated by ten cavity 
wavelengths or more. 

Finally, we address a crucial fabrication issue: Given that 



the best currently achievable deposition control for impuri- 
ties in Si is ^ 10 A |26|, the impact of small donor mis- 
placements on the proposed device operation must be ana- 
lyzed. A deviation 5r in the position of a donor from a max- 
imum of the cavity field (see Fig. [TJ introduces a variation 
Ari(7 ~ 2tt'^ {\5r\ / X)'^Q,c on the cavity vacuum Rabi fre- 
quency at the position of the donor. Here A = 7.8 /im 
is the cavity wavelength. This implies that \5r\ ~ lOOA leads 
to Afic 3 X 10 ~''^r^c'. The time for a typical two-qubit gate 
such as a/swap is t V^ij' leading to an error probability 
in this operation of p {AQ.c /^cY ~ 1 x 10^^, which 
means that our operation scheme is quite insensitive to rela- 
tively large (several lattice parameters) donor misplacement 
within the simple donor linear array architecture. This is in 
contrast with Kane's original exchange-based proposal, which 
leads to much more stringent conditions on impurity position- 
ing, and requires elaborate two-dimensional architectures to 
compensate for donor misplacement f2n\. 

Our estimations indicate that the present proposal could 
meet the conditions for a robust quantum computation device. 
The precise quantum control of atoms and ions in optical cav- 
ities, akeady demonstrated in several labs, and the fact that Si 
is the leading material in terms of processing and device fab- 
rication, with sophisticated techniques for impurity implanta- 
tion and high-Q microcavity construction, indicate that this 
system might be a viable candidate for a working quantum 
computer using achievable technological resources. 
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